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SUMMARY 


A supersonic triplet singularity has been developed which eliminates 
internal waves generated by panels having supersonic edges. The triplet is a 
linear combination of source and vortex distributions which gives directional 
properties to the perturbation flow field surrounding the panel. The theore- 
tical development of the triplet singularity is described together with its 
application to the calculation of surface pressures on wings and bodies. 
Examples are presented comparing the results of the new method with other 
supersonic methods and with experimental data. 


INTRODUCTION 


Significant advances have been made in the application of surface 
singularity techniques to the analysis and design of complex aircraft 
configurations in both subsonic and supersonic flows, (Refs. 1 and 2). The 
success of these methods is primarily due to the use of high-order source, 
vortex, or doublet singularities on surface panels, and by improvements in 
the formulation of the boundary conditions. 

The use of surface panel si ngul arities in supersonic flow introduces 
jumps in the perturbation velocity flow field across the Mach waves originat- 
ing along the panel leading and trailing edges. If these waves propagate 
into the interior of the wing or body, they can induce spurious perturbations 
on the opposite side of the configuration and seriously affect the strength 
of the surface singularity distribution. 

The supersonic triplet singularity is designed to eliminate the spurious 
internal waves generated by panels having supersonic leading edges. In par- 
ticular, the perturbation velocities in the "two-dimensional" region associ- 
ated with supersonic panel edges are exactly cancelled in the interior flow. 
This results in a well conditioned system of boundary condition equations and 
a correspondi ng improvement in the pressure distribution calculated on the 
exterior of the configuration. 

The application of the supersonic triplet concept to the aerodynamic 
analysis of bodies having arbitrary cross-section is reported in Ref. 3. 

This report describes the extension of the method to the analysis of swept, 
tapered wings, and includes details on the derivation of the wing and body 
triplet singularities. In addition, examples are presented comparing the 
results of this method with other supersonic methods and experimental data. 



LIST OF SYMBOLS 


a Panel inclination = tan 6 

Panel taper ratio 

b Panel span 

c Panel chord 

C Constant 

d Hyperbolic distance to panel corner 

0 Hyperbolic distance to panel apex 

F,G,H Velocity distribution functions 

£ Distance from panel apex to origin 

m Panel edge slope = tan 5 

M Mach number 

r,0 Angular coordinates 

u,v,w Perturbation velocities 

x,y,z Cartesian coordinates 

a Angle of attack 

3 Prandtl-Glauert factor =“\/m^ - 1 

Y Geometric parameter 

6 Panel inclination angle 

A Leading-edge sweep = tan A 

A Leading-edge sweep angle 

T 

Integration variables 


2 



Subscripts 


A Axi al 

C Ci rcumferenti al , 

E Edge 

L Line 

R Radial 

S Source 

T Triplet 


constant 


V 


Vortex 



AERODYNAMIC THEORY 


The triplet singularity is a linear combination of source and vortex 
singularity distributions. Since the vortex (or doublet) sheet may be con- 
sidered the result of combining, in a special limiting process, two source 
sheets of equal and opposite strength, the combination of a vortex sheet 

of unit strength with a source sheet of strength t = V , has been 

termed a triplet. The procedure is only applicable to panels having super- 
sonic leading edges, corresponding to real values for t. A simple illustra- 
tion of the basic concept for two-dimensional flow is given in Figure 1. In 
this example, the axial and normal velocity component vectors add in the flow 
field above the panel, but cancel exactly below the panel, resulting in the 
desired unidirectional perturbation velocity field. 

The extension of this concept to the analysis of three-dimensonal 
bodies and wings is described in the following sections. 


Body Triplet Singularities 

Three types of triplet singularities are required to apply this method 
to bodies having arbitrary cross-section. They are designated circumferential , 
edge, or radial, depending on the orientation of the elementary line singu- 
larities used in the derivation. An axial triplet singularity can also be 
obtained by combining the circumferential and edge triplets. 

Circumferential Triplet . The circumferential triplet is obtained by 
combining a circumferential source and vortex distribution. The elementary 
source and vortex filaments making up these singularities lie in the plane of 
the panel, parallel to the y-axis, as indicated on the following sketch. 

Note that the side edges of the panel intersect the x-axis a distance, £, 
ahead of the origin. 


X = 


y 



z) 


z) 
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TWO-DIMENSIONAL FLOW 



CONSTANT SOURCE CONSTANT VORTEX CONSTANT TRIPLET 

(LINEARLY VARYING (UN I -DIRECTIONAL) 
DOUBLET SHEET) 


cn 


Figure 1. - Supersonic Triplet Concept. 



The three components of velocity induced at a point, x,y,z, in the field 
by one corner of the panel are derived in Appendix A. The influence of the 
vortex distribution is given by Eqns. (A12) through (A14), and the source dis- 
tribution by Eqns. (A15) through (A17). The circumferential triplet is ob- 
tained by combining these singularities as follows. 


“'tc "" ^VC ^ ^SC 
= F + H - 3mG 


(1) 


'''tc ''vc * ®''sc 


= 6G 


(2) 


”'tc ^ ”vc ®“sc 

= 6(6mG - F - H) 

= -Su^P (3) 


where 


F = tan"'' - 3^2^ 

x(m(x + ii) - y) - B^mz^ 


(4) 


G = 


V7T 


cosh 


2m2 


3 m 


-1 X - 3^m(y - m£) 

3V(y - m(x + l))'^ + (1 - 3^m^)z^ 


(5) 


H = sin-’ -- 

Vx^ 


( 6 ) 
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Within the two-dimensional flow region (ahead of the Mach cones from the 
inboard and outboard corners of the panel, and behind the Mach wave from the 
unswept leading edge), F = tt above the plane of the panel, and -tt below, 
while H = n within this region on both sides of the panel. Thus the sum 
F + H appearing in Eqn. (1) and (3) has a value of 2tt above the panel, and 
zero below. This special relationship provides the desired wave cancellation 
property of the triplet singularity. 

The three components of velocity induced by a panel inclined at an angle, 

6, to the x-y plane are obtained by applying a Lorentz transformation about the 
y-axis. The geometry is illustrated on the following sketch. 


z Mach Wave 



For a panel lying in the x',y' plane, and defining a = tan 6 , 


u 


u ' - aw ' 

Vi - 


(7) 


V = v' 


( 8 ) 


_ w' - e^au' 

W = — i 

VI - 


(9) 


Expressions for the velocity components in the primed system of coordi- 
nates are given by Egns. (1) through (3). Substituting 
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(10) 


yi = X - e^az 

Vi - 
y' = y 

7* = z - ax 

V 1 - 

ml - m 

rn 

VI “ 


01 ) 

( 12 ) 

(13) 


The transformed velocity components in the reference coordinate system 
become: 


F + H - SmG 
^TC ' 1 - 3a 


V 


TC 


= G 


"tC = -S'^TC 


(14) 

(15) 

(16) 


with 


P = tan"^ (z - ax)Vx^ - B^(v - ml)^ - 

x[m(x + t) + y] + 6^2 [ay - m(z + at.)] 


(17) 
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G = 


1 

Vi - + m2) 


cosh"^ . — ~ - nr£) + az] 

3V[y^- m(x + £)]^ - 3^[ay - m(z + a£)]^ + (z - ax)^ 

(18) 


and H is given by Eqn. (6). 


Edge Triplet . The edge triplet is obtained by combining edge source 
and vortex distributions. The elementary source and vortex filaments making 
up these singularities lie in the plane of the panel parallel to one edge, 
as indicated on the following sketch. 



The three components of velocity induced at a point, x,y,z, in the field 
by one corner of the panel are derived in Appendix A. The influence of the 
vortex distribution is given by Eqns. (A24) - (A26), and the source distribu- 
tion by (A27) - (A29). In addition, the influence of a linearly varying line 
vortex along the leading edge (Eqns. (A57) - (A59)) is added to the vortex 
distribution, and a linearly varying line source along the leading edge 
(Eqns. (A63) - (A65)) is added to the source distribution. Combining these 
four contributions, the three components of velocity for the edge triplet 
singularity may be written: 


9 



u 


'tE "" “VE “VL * ®^“SE * “SL^ 

= I H - 6mG)] (19) 

v'te = ''VE ^ V + ®(''SE ^ ''sl) 

= F + H - BmG (20) 

”'tE " '*'VE ■*' ”VL ■*' ®^”SE ■'■ ”SL^ 

= G + Gm(F + H - 3mG) - ^ 

= G - (21) 

where F, G, an(J H are given by Eqns. (4) - (6), and 
D = Vx^ - 3^(y - mil)^ - 

Applying a Lorentz transformation by substituting Eqns. (19) - (21) in 
Eqns. (7) - (9), the three components of velocity induced at a point, x,y,z, 
in the reference coordinate system by an inclined panel can be written: 

“te " 6 [ X - sz ■ 1 -"'ea ''tE ■ 
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v-p^ = F + H - 6mG 


(24) 


w. 


TE 


= G - 3u 


TE 


(25) 


where F and G are given by Eqns. (17) and (18), H by Eqn. (6), and D by Eqn. 

( 22 ). 


Axial Triplet . The influence of an axial triplet distribution on an 
inclined panel can be obtained by adding the influence of an edge triplet 
to the influence of a circumferential triplet multiplied by the factor, Bm. 
Combining Eqns. (23) - (25) with Eqns. (14) - (16), the three components of 
velocity can be written: 


^TA ^TE ^TC 


.ir D 

e Lx - Bz 



(26) 


"tA = ''TE ^ S'" ''tC 


= F + H 


(27) 


"tA = WtE + '"TC 


= G - 



(28) 


This corresponds to a triplet distribution having the elementary source 
and vortex filaments parallel to the x-axis. 
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Radial Triplet . The elementary source and vortex filaments making up the 
radial triplet distribution lie inthe plane of the panel along radial lines 
from the intersection of the two side edges extended. This point lies on 
the x-axis a distance, £, ahead of the origin, as indicated on the following 
sketch. 


X = -£ 



£) 


X 




The three components of velocity induced at a point, x,y,z, by one corner 
of the panel are derived in Appendix A. The influence of the vortex distri- 
bution is given by Eqns. (A39) - (A41), and the source distribution by Eqns. 
(A57) - (A59). 

In addition, the influence of a linearly varying line vortex along the 
leading edge (Eqns. (A57) - (A59)) is added to the vortex distribution, and a 
linearly varying line source along the leading edge (Eqns. (A63) - (A65)) is 
added to the source distribution. Combining these four contributions, the 
three components of velocity for the radial triplet singularity may be 
written: 


“'tR ■ “VR “VL ®^“SR 

^ 3(x + Jl - Bz) [ X + 3z - X - 6z ®^^r] 

(29) 
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I 


'''tR ~ ''VR ''VL ^^''SR 


"" X + £ - 3z ^ " ^^^R^ 


( 30 ) 


'^'tR ^ '^VR ■*■ '^VL '^SL^ 

" X + £ - 6z [ X + "3z ^ ~ x“-~3z 

(31) 


where 


F IS given by Eqn. (4) , 

H is given by Eqn. (6), 

D is given by Eqn. (22), 

= j- rnch"^ x(x + £) - 3^v(y - m£) - 3^Z^ 


'R D 


I ■ •■■ - _ - — - _ 

R 6W (y - m(x + ^))^ + (1 - 3^m^)z^ 


(32) 


and 


Dp = V(x + ^)^ - 3^y^ - 3^z^ 


(33) 


Applying a Lorentz transformation, by substituting Eqns. (29) - (31) in 
Eqns. (7) - (9), the three components of velocity at a point, x,y,z, in the 
reference coordinate system by an inclined panel become: 


^TR " 3^ X + J - f{z + a£) ^ 


- ^yV-7 


- 3z 


+ 3(z + a5.)G 


R 


(34) 
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VjR = c(F + H - 3yGp) 


(35) 


w 


TR 


" ^ } X + £ - ^3'(z~+ a£) 




(36) 


, n - ga)Jl 
X + £ - 6(2 + a£) 

F IS given by Eqn. (17) 

H is given by Eqn. (6) 

D is given by Eqn. (22) 

Dr is given by Eqn. (33) 

G = J_ cosh~^ x(x + £) - 3^.v(y - m£) - 6^z(z + a£) 

^ ^R 6£V(y- rn(x + £))^ - 3^ (ay - m”(z^+ a£7)^ + (z- ax)^ 

(37) 

Rectangular Triplet . For the special case of a rectangular panel, 
m = 0 and £ -> oo. Taking this limit, 


where 


and 


and 


“tR ■ '^TE ' “TA 


1 

6 



(38) 


''tr ‘ ''te ‘ ''ta 


F + H 


(39) 
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(40) 


= G - 


D 


X - 3z 


The influence of an edge panel of a triplet panel group is obtained by 
subtracting the contribution of the radial triplet from that of the edge 
triplet. This has the effect of distributing the vorticity from a side 
edge to a trailing edge of a panel. For a rectangular panel, a special 
limiting process is required to obtain the correct influence. 


Uy = lim 


[jI(Ut£ - Uyj^)j 


“ y(F + H) + zG + D/3 


(41) 




lim 

£ ^ 00 ^ 


TR 


'] 


= (x - 3z)(F + H) + 3yG 


(42) 


w. 



w 


TR 


>] 


= (x - ez)G - 6y(F + H) - 2D 


(43) 
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Wing Triplet Singularities 


In the analysis of wings, the triplet singularity is required only for 
those panels having supersonic edges. For panels having sonic or subsonic 
edges, a combination of constant sources and linearly varying vortices, as 
described in References 2 and 4, provides a suitable aerodynamic representa- 
tion. The influence of a constant triplet singularity, located on a panel 
having a swept supersonic leading edge is derived as follows. 

Constant Triplet . The constant triplet is obtained by combining a 
constant source and vortex distribution. The elementary source and vortex 
filaments making up these singularities lie inthe plane of the panel parallel 
to the leading edge, as indicated below. 



The three components of velocity induced at a point, x,y,z, in the field 
by one corner of the panel are derived in Appendix B. The influence of a 
constant vortex distribution is given by Eqns. (B7) - (B9), and a constant 
source distribution by (B13) - (B15). The constant triplet is obtained 
by combining these singularities as follows. 


^TC ^VC IJ5C 

= F + H - m Vb" - G 


(44) 
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I 


''tc = ''vc ^ Vsc 


= -X(F + H) + G 


(45) 


w 


TC 


= w 


VC 


V 


3^ - w 


SC 


= -Vs^ - (F + H) - (X - 6^m) G 


(46) 


where F and G are given by Eqns. (BIO) and (BIT), and 


H = cos 


Xx - B^y 

bV(x - Xy)^ + (X^ - B^z^) 


(47) 


The influence of the complete panel is obtained by combining the in- 
fluences of each of the four corners. Within the "two-dimensional" region 
associated with the supersonic leading edge, the sum of F and H is 2 it above 
the plane of the panel, and zero below. This special relationship again 
provides for cancellation of the perturbation velocities in the interior 
of the wi ng. 

A pair of linearly varying line vortices is added along the side edges 
of the panel to transfer the accumulated vorticity into the wake. The 
influence of these line vortices is given by Eqns. (Ill) - (113) of 
Reference 2. 
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APPLICATION 


No special procedures are required to apply the triplet singularities 
to the analysis of wings and bodies in supersonic flow. One control point 
is associated with each panel (or group of panels) used in the aerodynamic 
representation. The boundary condition of tangential flow is then imposed 
at each control point, and the resulting system of linear equations is 
solved to determine the individual singularity strengths. Surface pressures, 
forces, and moments acting on the configuration are then calculated in the 
usual way. A computer code has been developed to perform the calculations 
and is available as a modification to the USSAERO program. 


Body Panel Arrangement 

A panel arrangement suitable for the analysis of a body having arbitrary 
cross-section is shown in Figure 2. The panel leading and trailing edges are 
defined by planes perpendicular to the reference x-axis, while panel side 
edges are defined by the body meridian lines. A grouping of six panels is 
used to make up a triplet singularity which will satisfy the Helmholtz vortex 
conservation laws without requiring line vortices along panel edges or trail- 
ing vortex wakes. Within this six-panel group, the individual singularity 
strengths are prescribed to ensure that the elementary vortex filaments form 
closed loops, as indicated on the figure. The center row of panels contains 
only circumferential triplet singularities of equal and opposite strengths, 
ihe singularities in the outer rows are composed of a special combination of 
edge and radial triplets, which result in zero vorticity along the outer two 
edges, and match the inflow (or outflow) of vorticity along the inner two 
edges. 

The strengths of the triplet singularities in each group are determined 
by satisfying Neumann boundary conditions at the control points. 


Wing Panel Arrangement 

A conventional panel subdivision is used on the wing. The panel side 
edges are defined by planes parallel to the reference x-axis, while the 
leading and trailing edges are defined by constant percent chord lines. 
Panels are located on both upper and lower surfaces of the wing. 

For wings having subsonic leading and trailing edges, a combination of 
constant sources and linearly varying vortices are used in the aerodynamic 
representation. The linearly varying vortex singularities are distributed 
over two adjacent panels, as illustrated in Figure 3, and include a pair of 
constant vortices trailing downstream along the panel edges and into the 
wake. The principle of symmetrical singularities is applied, as described 
in Reference 5, by equating the source and vortex strengths on corresponding 
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SECTION A-A 


Figure 2. - Panel Arrangement on Arbitrary Body. 




I 


panels on the upper and lower surfaces of each section. The details of the 
extension of this method to the analysis of three-dimensional wings is 
given in Reference 4. 

For wings having supersonic leading and trailing edges, constant triplets 
alone are used on each panel. This aerodynamic representation minimizes inter 
ference effects between the upper and lower surfaces of the wing. 

Wings having subsonic leading edges and supersonic trailing edges present 
a more difficult problem. In this case, all panels ahead of the sonic genera- 
tor line have subsonic edges, and those behind this line have supersonic 
edges. (A sonic generator line is defined as a wing generator line parallel 
to the Mach line.) A combination of constant sources and linearly varying 
vortex singularities are used ahead of this line, changing to constant triplet 
distributions on panels behind the line, with the linear distributions over- 
lapping the constant distributions on the panel containing the sonic genera- 
tor line itself. The strengths of the singularities are then determined by 
applying the symmetrical singularity method described above. This technique 
results in relatively smooth variations in singularity strength across the 
sonic generator line. 
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EXAMPLES 


Several examples are presented to illustrate the application of the 
triplet singularity method to a variety of wings and bodies in supersonic 
flow. Theoretical pressure coefficients are calculated using the isentropic 
pressure coefficient formula. The examples are compared with other theoreti- 
cal methods and experimental data. 


Parabolic Body of Revolution 

Figures 4 and 5 give the axial pressure distribution on a parabolic body 
of revolution having a fineness ratio of 20 at M = /2, Comparisons between 
the triplet method and the surface source method of Reference 2 are shown 
for ct = 0 and 5°. For this smooth slender body, the two methods give very 
similar results, although a close examination shows small irregularities in 
the pressure distribution calculated by the source method caused by internal 
wave reflections. 


Cone-Cyl i nder-Cone 

A more severe test case was provided by a 15° cone-cylinder-code body at 
M = 2.0. In this example, the surface source panel method of Reference 2 
fails to give a convergent solution. Figures 6 and 7 show the axial pressure 
distribution at oi = 0 and 5°, respectively, compared to that obtained by the 
classical Karman-Moore method. (The Karman-Moore method uses a distribution 
of line sources and doublets along the axis, together with tangency boundary 
conditions on the body surface. Ref. 6.) The results of these two methods 
agree closely, indicating that the triplet method was effectively eliminated 
the strong internal waves that caused the source method to diverge. Figure 
8 shows the circumferential pressure distribution on the 15° cone at M = 2.0 
and a = 10°, on a section just behind the shoulder of the cone-cylinder and 
on a section of the cylinder extended 35 diameters behind the shoulder. The 
pressure distribution on this last section is seen to approach the incompres- 
sible cross-flow on a circular cylinder in two dimensons. 


Nacelle with Internal Flow 

The exterior pressure distribution on a circular nacelle with internal 
flow has been calculated at M = for a = 0 and 10°. The results are pre- 
sented on Figure 9. This example again indicates the effectiveness of the 
triplet singularity in suppressing spurious internal waves, and can be com- 
pared with the nacelle results presented in Reference 1. 
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Figure 5. Parabolic Body; n = 5 
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Figure Circumferential Pressure on 15° Cone-Cylinder. 
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Figure 9. Nacelle Pressure Distribution 






Ogi ve-Cyl i nder-Boattai 1 

Figures 10 and 11 give the axial pressure distribution on an ogive- 
cylinder-boattail body at M = 2.3. For a = 0°, Figure 10, a comparison is 
made between the triplet method and the experimental data of Reference 7. 

The theory tends to underestimate the pressure on the nose of the body, but 
agrees well with the data elsewhere, except on the boattail where viscous 
effects dominate. The underestimation of the nose pressure is a result of 
applying the linearized potential flow equations and tangency boundary con- 
ditions. A better approximation to the pressure on the nose can be obtained 
by applying a local Mach number correction described in Reference 4. The 
improvement is indicated by the dashed line on the figure. 

For a - 4°, Figure 11, a comparison is made among the triplet method, 
the finite-difference method. Reference 8, and experimental data. The 
finite-difference method agrees well with the experimental data except on 
the boattail, while the triplet method again underestimates the pressure on 
the nose. No local Mach number correction has been added to the triplet 
results in this case. 


Elliptic Cone 

The triplet singularity may also be applied to the analysis of non- 
circular bodies. In Figure 12, the circumferential pressure distribution 
on an elliptic cone at M ~ 1,89 is compared with experimental data. Reference 
9. For a = 0°, the theory overestimates the pressure slightly. For a = 6° , 
a more pronounced overestimation occurs on the lower surface. 


B-1 Forebody 

The pressure distributions along the upper and lower meridians of the 
B-1 forebody for M - 2.2 and a = 2° are shown on Figure 13. Results from the 
triplet method are compared with the finite-difference method and experi- 
mental data obtained from Ref. 10. The finite-difference method agrees 
exceptionally well with the experimental data in this example. On the other 
hand, the triplet method underestimates the pressure bn the nose cone, and 
fails to predict the pressure peak behind the canopy shock. This result is 
due to the shortcomings of linearized potential flow theory and can only be 
partially offset by applying a local Mach number correction. 


Uncambered Arrow Wing 

The pressure distribution calculated for an uncambered arrow wing having 
70° sweepback and a 3 percent biconvex section at M = 2.05 and a = 4° is 
presented on Figure 14. The results are compared with experimental data from 
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Fiqure 13. B-1 Forebody. 





Figure 14. Uncambered 


Arrow Wing; M - 2.05, a = 4 Degrees. 
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Reference IT, and a previous USSAERO calculation using planar boundary condi- 
tions with linearized pressure coefficient formula. 

This wing has a subsonic leading edge and a supersonic trailing edge, 
and consequently uses a combination of linear vorticity and constant sources 
ahead of the sonic generator line, and constant triplets aft of that line. 

The results are smooth and are generally in good agreement with experimental 
data except on the upper surface near the wing tip. The theory also fails to 
predict the development of a small vortex above the leading edge on the in- 
board sections. The nonplanar boundary condition applied in this example 
gives a noticeable improvement in the pressure distribution on the lower sur- 
face, compared with the previous planar boundary condition calculation. 
However, much of this improvement is due to the use of the full isentropic 
pressure coefficient formula rather than to the particular aerodynamic model 
empl oyed . 


Cambered Arrow Wing 

Figure 15 gives the pressure distribution calculated for a cambered and 
twisted arrow wing having the same planform and airfoil section as the pre- 
vious example. The camber and twist was selected to give a design lift 
coefficient of .08 at M - 2.05 (see Ref. 11). The theory is compared with 
experimental data for a = 4°, and also with a previous USSAERO calculation 
using the planar boundary condition option. 

In this calculation, the effect of the camber and twist was obtained by 
analyzing the equivalent uncambered wing in a nonuniform onset flow consist- 
ing of the normal component of the free stream velocity, plus an incremental 
normal velocity proportional to the slope of the chordwise camber distribu- 
tion. The boundary condition of tangential flow is then applied at the 
upper and lower surface control points on the equivalent uncambered wing. 

This is consistent with the technique used to add camber effects in the planar 
boundary condition option of the USSAERO program. 

The resulting chordwise pressure distributions are in reasonable agree- 
ment with the experimental data, except on the lower surface near the wing 
tip. Again, the nonplanar boundary condition calculations agree more closely 
with the lower surface experimental data than the previous planar boundary 
conditon calculations, due primarily to the pressure coefficient formula 
employed. 
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Figure 15. Cambered Arrow Wing; M = 2.05, a = 4 Degrees. 
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Figure 15. Concluded. 


CONCLUSIONS 


A triplet singularity has been developed and applied to the analysis 
of wings and bodies in supersonic flow. The new singularity has been found 
most effective in the analysis of bodies, and appears to eliminate any 
sensitivity of the surface pressure to the details of the panelling. 

The new singularity also makes it possible to extend the method to the 

analysis of wings defined by surface panels, using nonplanar boundary condi- 
tions in supersonic flow. However, no significant improvement in the predic- 

tion of the wing pressure distribution has resulted from this extension. It 
is concluded that a more extensive evaluation of the method is required to 
determine if the overall improvements obtained justify the additional numeri- 
cal complexity involved. 
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APPENDIX A: DERIVATION OF BODY SINGULARITIES 


1. Circumferential Vortex Distribution 


The incompressible velocity components induced at a point, P, by a constant 
vortex distribution in the plane of the panel are derived by integrating the 
influences of a series of elementary line vortices extending across the panel 
parallel to the leading edge. The geometry of an elementary line vortex 
located at a distance, C, from the leading edge and having strength, d^, is 
illustrated below. 



The panel lies in the x,y plane and the edge slope, m = tan 0, is an 
arbitrary constant. The distance from the point, P, to the left end of the 
elementary vortex is 


d = V (x - + (y ” nix)^ + 


The velocity components are obtained by applying a 90 degree coordinate 
rotation to the line vortex velocity formulas given by equations (6) through 
(8) of Reference 2, and integrating across the panel as follows. 
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(AT) 



/ 


(y - mg) dg 

({x - 5)^ + z^)V(x - ?)^ + (y - mzV + 


= 0 


(A2) 


w. 


/ 


(x - c)(y - m^) d? 


((x - + z2) V(x - + (y - m ^)2 + 


(A3) 


Put V = X - Then 


= z [(y - mx) + mJ^] (A4) 

= mz^ - (y - mx) J 2 - mK^ (A5) 


where 



d v 

+ z^) “\/av^ + 2bv + c 


V dv 

+ z^) + 2bv + c 

dv 

^ + 2bv + c 


(A6) 


(A7) 


(A8) 


The integrals, J*| and J 2 » are given in Appendix I of Reference 2, and the 
integral K-j is in standard form. Evaluating these integrals, the three 
components of velocity induced by the inboard corner of the panel can be 
wri tten : 
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I 


(A9) 


-1 zVx^ + 

u = tan — 

-x(y - mx) + mz^ 



(A10) 


w ^cinh-1 — y_.. x + my 

c yx^ + z^ yl + y(y - mx)^ + (1 + m^) 


(All) 


The compressible velocity components can be obtained from the above 
expressions by applying the extended Gothert's rule. In supersonic flow, with 

3 = 1, 


, zVx^ - 3My^ + Z^) 

u = F = tan — (A12) 

c -x(,y - mx) - m3^z^ 


V 


V 


c 


= 0 


(A13) 


w„ = B 


cosh-^ 


Vi - 


X - B^my 


B V(y - nix)2 + (1 - B^m^) 


- sin 




Vx - 


= B fern G - H) 


(A14) 
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2. Circumferential Source Distribution 


The incompressible velocity components induced at a point by a constant 
source distribution in the plane of the panel are given by Eqns. (18) through 
(20) of Reference 2. In supersonic flow, these expressions become: 

= - (3m G - H)/3 (A15) 

V3, = G (A16) 

W3, = -F (A17) 

where F, G, and H are defined above. 



The incompressible velocity components induced at a point by a constant 
vortex distribution oriented parallel to the side edge of the panel are derived 
by integrating the influence of elementary line vortices across the panel in 
the spanwise direction. The geometry is illustrated below. 
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The panel lies in the x,y plane, and the edge slope, m= tan e, is an arbi- 
trary constant. The distance from the point, P, to the origin of the vortex 
is 


d = + (y - ri)^ + 


The velocity components are obtained by rotating the line vortex 
velocity formulas given by equations (6) through (8) of Reference 2 through 
the angle, 6, and integrating across the panel as follows: 


^VE 



x + m(y - n) dn 

[(y - mx - n)^+ (l+m^)z^]Vx^ + (y - q)^ + 


(A18) 


Vy^ = - u/m 


(A19) 





(y - mx - T]){x + m(y - n)) dn 
[(y “ mx - n)^+ (l+m^)z^]Vx^ + (y - n)^ + z^ 


(A20) 


These integrals can be expressed in terms of the J and K integrals by sub- 
stituting V " y - mx - n. Performing the integrations, the three components 
of velocity induced by the inboard corner of the panel may be written: 


^VE = 


-m tan 


zV x^ + y^ + z^ 

-x(y - mx) + mz^ 


(A21) 


V 


VE ' 


-u/m 


(A22) 
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= V 1 + sinh~^ - > -. - -.- A , ^y . 


V (y - nix)2 + 0 + m^)z2 


m sinh ^ ~=^== + m log 

V + 2^ 


2 + 2 ^ 


(A23) 


These formulas can also be obtained by performing a 90° coordinate rotation in 
Eqns. (43) through (45) of Reference 2. , In supersonic flow. 


'VE 

(A24) 

VE = ' 

(A25) 

“ (1 - 3^m^)G + BmH 

(A26) 


where F, G, and H are defined by Eqns. (A12) through (A14). 


Edge Source Distribution 

The corresponding edge source distribution can be obtained by multiplying 
the ci rcumferential source distribution by -m. In supersonic flow. 


u^^ = m(BmG - H)/B 

(A27) 

''SE " 

(A28) 

“SE ” 

(A29) 
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5. Radial Vortex Distribution 


In the radial vortex distribution, the velocity components at a point are 
obtained by integrating the influences of a series of elementary line vor- 
tices lying along radial lines from an apex defined by the intersection of the 
two side edges. The geometry is illustrated below. 


x = 


// 


/ 

/ / 
/ / 

/ / 
f / 


n 

/ 


-I 



The elementary line vortex lies along the line, y = mx, and originates 
at the point, o, mi?-, where I is the distance of the apex ahead to the panel 
leading edge. The distance from the point, P, to the origin of the vortex is 

D = V + (y - m£)^ + 


In this case, the vortex inclination, m = tan 0 , varies from m-| to m 2 
across the oanel. The velocitv comoonents are obtained bv rotating the line 
vortex velocity formulas given by Eqns. (6) through (8) of Reference 2 through 
the angle, 9, and integrating with respect to m as follows: 
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m. 


^VR ~ 


/ 


[x + m(v - dm 




[(y - ni(x + Z))^ + (1 + m^)z^]Vx^ + (y ~ mzy + 


(A30) 




VvR = -zZ 


/ 


X + m(y - m^) dm 


m. 


[(y - m(x + Z))^ + (1 + + (y - ni5,)^ + z^ 


(A31) 


^vR + yv] /z 


(A32) 


vfx + 5,^ 

Substituting v = m + — , the above integrals can be expanded 

(x + Z)^ + z^ 

into expressions containing the J integrals (with y = yz£/((x + £)^ + z^), 
and other integrals of standard form. Performing the integrations, the 
three components of velocity induced by the inboard corner of the panel 
may be written: 


'VR 


-1 S 

(x + SL)^ + I 


[ D + £(x + J.)Gpl+ - zf p 

^ J (x + Z)^ + z^ 


2yziix_i.y_ [yG 

(x + ly + L R 



(A33) 
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''VR 


£ 

(x + £)^ + 2^ 


(x + £) F + z 



(A34) 


Wwp = < (x + £)(D + £(x + £)Gp) 

(x + £)2 + zM ^ 


+ yj t(x + - . f (yg . H) . 2yz^( . x + »■) , F (A35) 

(x + £)^ + ^ (x + £)^ + ) 


where 


F = tan 


zD 

-x(y “ m(x + £)) +mz^ 


(A36) 




i sinh'^ x(x + g.) + .y(y - ttiii-) + 

t V{y - ni(x + l))^ + (1 + m^)z^ 


(A37) 


H “ s1nh 


V - tn£ 

V 


(A38) 


and 


D = Vx^ + (y - fn£)^ + z^, 


Dp = V (x + 
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In supersonic flow, 


where 


-1 


‘VR 


(x + £) - 


j z(D + £(x + J)Gp)+ yg.-^ + F 

( ^ (x + g)^ - 


_ 2 Syzg(x g) ( ye . H) 
(x + g )2 - g^z^ 


(A39) 


r„p = ^ I (x + g)F - 6 z(ByGp - H) ( 

VK q2^2 ' 


(x + i)^ - B^z 


(A40) 


w 


1 


VR 


(x + - B^z^ 


(x + £) (D + £(x + ^)Gj^) 


-m 


(x + + B 


2z2 


2B^VZil(x + £) 


(ByGn - H) + ' - • / — F 

(x + £)" - B"z^ ^ (x + £)" - B^z" 


(A41) 


F = tan 


zD 

-x(y - m(x + il)) - B^mz^ 


(A 42 ) 


^R 


±_ x(x + l) - B^y(y - mil) - B^z^ 

^R B^V(y - tn(x + £))^ + (1 - B^ni^)z2 


(A43) 
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6. Radial Source Distribution 

The derivation of the velocity components induced by a radial source 
distribution is similar to that used for the radial vortex distribution, 
with elementary line sources replacing the elementary line vortices along 
radial lines from the apex of the panel. The three components of velocity 
may be written in integral form as follows: 
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V 


SR 



[ xy - m(x(x + £) + z^)] mdm 

- m(x + £))^ + (1 “ m^)z^]Vx^ + (y - m£)^ + 


(A46) 


'^SR " '^R 


(A47) 


Carrying out the integrations, using the substitutions indicated in the 
previous section, the three components of velocity induced by the inboard 
corner of the panel in incompressible flow are: 


‘SR 


1 

(x + ly + T 


(x + £)D + zHGp - — 

^ (x + ly + z^ 


+ 


(x + ly - z^ 
(x + ly + 




(A48) 


zF - (x + Sl)(yG - H) - H (A49) 

(x + Jl)2 + zM ) 


'^SR “^VR 


(ABO) 


where F, Gr, and H are defined by Eqns. (A36) through (A38). 
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In supersonic flow; 


S^[(x + «)" - 6 "z "3 


(x + «.)D - - 


2e^vzi(x + );,) p 
(x + i)^ - 


+ gyx (x + zy + e^z=^ 


(3yGp 



{A51) 


'SR 


3[(x + £)^ - 


j I 6zF - (x + )l)(6yGp - H) I - ^ 


(A52) 


“SR ■ ~'^VR 


(A53) 


where Um is defined by Eqn. (A39), and F, Go and H are defined by Eons. (A42) 
through (A44). 


7. Leading-Edge Line Vortex 

In incompressible flow, the three components of velocity induced at a 
point, (x,y,z), by a linearly varying line vortex located along the panel 
leading edge are: 


zd 


'VL 


X^ + Z' 


0 


VvL = 0 


(A54) 

(A55) 
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w 


VL 


xd 

+ 2 ^ 


0 


where 


d = 


In supersonic flow. 


zd 




VvL = 0 


(A56) 


(A57) 


(A58) 


where 



d = \/^ 


2<d 

- 


- 3^(y^ + 2^) 


(A59) 


8. Leading-Edge Line Source 

In incompressible flow, the three components of velocity induced at a point, 
(x,y,z) by a linearly varying line source located along the panel leading edge 
are : 


u 


SL 



(A60) 


"SL 


si nh 



(A61) 
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where 



In supersonic flow, 

U = xd 

- e^z^) 



(A62) 


(A63) 


(A64) 


(A65) 
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APPENDIX B: DERIVATION OF WING SINGULARITIES 


1 . Constant Vortex Distribution 

The velocity components induced at a point, x,y,z, by a constant vortex 
distribution on an unswept body panel in supersonic flow are given by equa- 
tions (A12) through (A14). The velocity components corresponding to a swept, 
tapered wing panel are obtained by applying a Lorentz transformation about the 
z-axis. If the leading edge of the panel is swept back through the angle. A, 
and X = tan A, the three components of velocity may be written: 


“vc “vc ■ 

(B1) 

''VC = ''VC - 

(B2) 

Wyc = V B" - x^ w^^/6 

(B3) 


where u' , 
VC 


v' , W' . 
VC VC 


are given by Eqns. (A12) through (A14), with the following 


coordinate transformations applied. 


= B(x - Xy) 

Vs^ - 


(B4) 


Sy' 


B^V - Xx 

Vb^ - X^ 


(B5) 


Bm' = 

B (1 - mX) 


(B6) 


where m = tan 0. 
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The geometry is illustrated in the following sketch: 



After some simplification, the velocity component formulas for a constant 
vortex distribution on a panel having arbitrary edge sweep in supersonic flow 
may be expressed as follows: 


'vc " '' 

(B7) 

VC " 

(B8) 

-y(, = - e^)H - {A - 6"m)G 

(B9) 


where 


F = tan"^ 0 - z ~ + z^) 

(x - Ay)(mx - y) + (A^ - 3^m)z^ 


(BIO) 


G = 


vr^ 


cosh 


-1 


2m2 


3 m 


X - 3^my 

3V(y - nix)^ + (1 - 3^ni2)z^ 


(BID 
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and 


1 


for X > 3 


H = 


V 


cosh 


-1 


- 3* 


Xx - 3^.v 

3V(x - Xy)^ + (X^ - B^z^) 



cos 


Xx - 3^V 

bV(x - Xy)^ + (X^ - 3^z^) 


for X < 3 


(B12) 


It should be noted that, for m = 0, Eqns. (B7) through (B9) reduce to 
Eqns. (108) through (110) of Reference 2. 


2. Constant Source Distribution 


The velocity components induced at a point, x,y,z, by a constant source dis- 
tribution on an unswept, tapered body panel in supersonic flow are given by 
Eqns. (A15) through (A17). The corresponding formulas for a constant source 
distribution on a wing panel having arbitrary edge sweep in supersonic flow 
can also be obtained by applying the Lorentz transformation, as described 
above. The resulting expressions are: 


1 ^^ = H - mG 

(B13) 

SC = G - XH 

(B14) 

'sc 

(B15) 


where F, G, and H are defined by Eqns. (BIO) through (B12). For m = 0, the 
above velocity component formulas reduce to Eqns. (96) through (98) of 
Reference 2. 
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3, Constant Edge Vortex Distribution 


The constant edge vortex distribution is obtained by combining the influ- 
ence of a radial vortex distribution on a swept, tapered panel with the 
constant vortex distribution derived for the same panel. Expressions for 
the three components of velocity, u^^, v^^, w^^, induced by the constant vor- 
tex distribution are given oy Eqns. 008) through (110) of Reference 2. Ex- 
pressions for the velocity components induced by the radial vortex distribu- 
tion are obtained by integrating the influence of an elementary line vortex 
of strength, d^, located along a radial line passing through the intersection 
of the leading and trailing edges extended. The geometry is illustrated in 
the following sketch. 



P(x,y,z) 


In integral form: 


c 



0 


(B16) 

(B17) 

(B18) 
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where 


X(x - 5) + y 

V (x - 5)^ + y2 + z2 


= (x - ^ “ Ay)^ + (1 + A^)z^ 


and A = A^ - aC/c 


Here, a = A-| - A 2 , b = c/a, and A-j and A 2 are the slopes of the leading and 

trailing edges, respectively. Making the substitution given by Eqn. (51) of 
Reference 2, and integrating, the velocity components may be written as 
follows: 


'VR 


VR 


.1. ( 

p2 I 


(c - ay) F, - azG 


= .i ( 


(cA - ax) yF-, 

P^- i 

+ az [aD/c - (c - ay)(3^G^ + 2{cA - ax)G^/p^)] 


(B19) 


+ z(x + D)/cr^ 


(B20) 


w 


VR 


1 


(cA 


P^ ( 

- 2az(cA 


ax) - (c - ay) [aD/c - B^(c - ay)G^ 
ax)F-j/p2] I - y(x + D)/cr^ 


(B21) 
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where 


^ = (c - ay)^ - 
(c - ay)^ + 


D = V (x - 5)^ + 3^r^ 


r^ = y^ + z^ 


32 = 1 - M‘ 


- (cX - ax)G^ - G 2 


and the remainder of the terms are defined by Eqns. (118) through (120) of 
Reference 2, It should be noted that the influence of a linearly varying 
line vortex has been added along the side edge of the panel. 

The constant edge vortex distribution is obtained by combining the above 
expressions as follows: 


^VE ~ ^^E ^VR " 

VyE = Vyj^ - w^^)/a 


Wy^ (c^ - w^^)/a 


(B22) 

(B23) 

(B24) 


where c^ is the chord length of the edge. 

A special case is required for swept constant chord panels (a = 0). 
this case, 

^VE " '^’^1 “ ^^^^1 " ^ 2 ^ 

Vy^ = (x - 2Xy)F^ + 2 ( 3 ^ + 2A^)G^ 


w 


VE 


= 2AzF^ + y((3" + A2)G^ - AG 2 ) - (x - Ay)(AG^ - - 2D 


In 

(B25) 

(B26) 

(B27) 
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4. Constant Vortex Wake 


The velocity components induced by a constant vortex wake are obtained by 
integrating the influence of elementary line vortices located along lines 
parallel to the x-axis across the panel in a spanwise direction. In incompres- 
sible flow, 




(B28) 


Avyj. = z f 

= . f h . - n )dn 

™ J {y- n)^ + 


(1 M 3 .) 

0 


(B29) 


(B30) 


where d = V(x - ^ - Xn)^ + (y - n)^ + 


Performing the integration. 


^''VW " *"l " ^2 

(B31) 

AWvy ~ ~ XG-| + G 2 ~ log r 

(B32) 

where , F^, G^ and G 2 are defined by Eqns. (99) through 

(102) of Reference 2. 

In supersonic flow. 


AVvw = Ft 

(B33) 

^'^VW ~ ~ ^^1 ^ ^2 

(B34) 


where G;| and are now defined by Eqns. (103) and (104) of Reference 2 . 
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5. Linearly Varying Source Distribution on a Tapered Panel 

The three components of velocity induced at a point, P(x,y,z), by a linear- 
ly varying source distribution in the plane of a swept, tapered panel are 
derived by combining the influences of elementary constant and linearly 
varying line sources located along radial lines passing through the inter- 
section of the panel leading and trailing edges extended. The geometry is 
illustrated in the fol owing sketch: 



y 


U 


+ 


Xb,b) 


The velocity components are obtained by rotating the line source veloc- 
ity formulas given in Reference 1 through the angle, A, and integrating across 
the panel chord as follows. 

For a constant line source. 


' f 


VT 


f (x - g - X.y)(X(x - g) + y) 

J (1 + rH J 


X 


0 + X^)d 


- dg 


(B35) 
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(B36) 


''C = 




f 


dg 

(1 + X^)d 


/ 


A(x - g - Ay)(A(x - g) + y) 
(1 + A^) r^d 


dC 


Wp = Wp = z f ~ ^ 

^ ^ J p2j 


(B37) 


where A = tan A = A^ - a?/c 

a = A^ - A^ 

d = V (x - + y 2 + 2 ^ 

= (x - 5 - Ay)^ + (1 + A^)z^ 

These formulas may be compared with Eqns. (31) through (33) of 
Reference 2 . 

For a linearly varying line source. 



V 1 + 
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X 


f (x - g- Xy)d 

J {^ + X2)r^ 


d? 


■/ 


0 + X2)3/2 


sinh 


^(x - g) + V 


d? 


(B38) 


V 


L 


AVj^ + u 

Vi + X 


L 

2 



X(x - ^ - Xy)d 
0 + X^)r2 



1 

0 + 


sinh”^ 


j^X X - g) + y (jr 
r ^ 

(B39) 


w 


L 



(B40) 


The influences of the constant and linearly varying line sources are 
multiplied by ^ and combined in the following manner to obtain the final 
integral forms. 


“SL = “C - 


■/ 


(x - C - Xy) [(cX-| - ax)(x - C) + y(c - ay) - az^] 

(1 + X^)r^cd 


m 


/ttt 


. f + 1 / XL_^ sinh-'' - - ^ y d5 

J 0 + X=^)d (1 + 
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Noting that 


1), and simplifying. 


form. 


aXg 


1 /■' + 


(1 + \ 1 + 


“SL ^ c/ ^ 


sinh 


-1 X (x - g) + y 


d? 


^ s i n h ^ 


X 2 (x - c) + y 


VTTT7 


V(x - c - X,y)^ + (1 + x„^)z' 


(B41) 


The integral appearing in this expressions cannot be evaluated in closed 
It is integrated numerically in the computer program. 


''SL = \ 


- f- ^ 1 f 

^ n + x2)d (1 + x^) 


3/2 


,-nh-'' Mx- - g) + y 


sinh 


dC 


-/ 


X(x " ^ - Xy) [(cX^ - ax)(x - C) + y(c - ay) - az^] 

(1 + X^)r^cd 




1 ( / 2- -1 y 

= " { Vl + X, sinh / ■ • 

M 1 V(x - X^y)^ + (1 + X^ 


2)z2 1 2 
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(B42) 


[^(c - ay)(e2G^ - 


(cX-| - ax)G2 + aZ(cA^ 



WsL = W(, - aw^/c 


z 

c 



- ax)(x - ?) + y(c - ay) - az 



= — jsF, + z 


fe^G^ - (cX. 



(B43) 


In incompressible flow, , G^ , G^, s, and p are defined by Eqns, (56) 

through (58) of Reference 2. in supersonic flow, corresponding expressions 
for these functions are defined by Eqns. (118) through (120). 

For a = 0, the above velocity component expressions reduce to Eqns. (37) 
through (39) of Reference 2. 
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